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Abstract 

We consider thermal strongly-coupled Af = 2 SYM theory with fundamental 
matter at finite isospin chemical potential. Using gauge/gravity duality, i. e. a 
probe of two flavor D7-branes embedded in the AdS black hole background, we 
find a critical temperature at which the system undergoes a second order phase 
transition. The critical exponent of this transition is one half and coincides with 
the result from mean field theory. In the thermodynamically favored phase, a 
flavor current acquires a vev and breaks an Abelian symmetry spontaneously. 
This new phase shows signatures known from superconductivity, such as an 
infinite dc conductivity and a gap in the frequency-dependent conductivity. 
The gravity setup allows for an explicit identification of the degrees of freedom 
in the dual field theory, as well as for a dual string picture of the condensation 
process. 
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Gauge/gravity duality is a powerful tool to compute observables of field the- 
ories at strong coupling. The original duality for superconformal field theories 
has been extended to field theories at finite temperature [l[ , using an AdS black 
hole as gravity dual. Moreover, the original correspondence has been extended 
to field theories with flavor degrees of freedom in the fundamental representation 
of the gauge group and to finite densities. In 0, a finite 517(2) isospin density 
d was considered. At high isospin densities d > d c , this configuration becomes 
unstable against flavor current fluctuations, corresponding to vector mesons Q. 
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Figure 1: Phase diagram for fundamental matter in thermal strongly-coupled J\f = 2 SYM 
theory Q], with fi the isospin chemical potential, M q the bare quark mass, M = 2M q A~ 1 / 2 , 
A the 't Hooft coupling and T the temperature: In the blue shaded region, mesons are stable. 
In the white and green regions, the mesons melt. Here the new phase is stabilized while it 
was unstable in Q{. In this phase we find some features known from superfluids. 

The central result of the present letter is that there is a new stable phase 
characterized by a g-meson superfluid. We show explicitly that the new phase 
shown in fig. 1 is thermodynamically preferred. Moreover, we check that the 
phase is stable (at least against the fluctuations which caused the instability 
beforehand) . 

Moreover we find that this phase has properties of a superconductor, such 
as a second order phase transition with a critical exponent of 1/2. The fla- 
vor current, which is analogous to the electromagnetic current, displays infinite 
dc conductivity and a gap in the frequency-dependent conductivity. Further- 
more, we give a dual string-theoretical picture of the Cooper pairs and find a 
dynamical generation of meson masses similar to the Higgs mechanism. The 
stabilization mechanism is motivated by the results of 0, [E[ ■ In the context of a 
phenomenological AdS/QCD model, a similar stabilization mechanism was first 
used in [6| to obtain a p-wave superconductor. 

Our results provide an important link between gauge/gravity with flavor and 
the recent gauge/gravity realizations of superconducting states e.g. in [1, 0, S|, 
as well as a realization of a superconducting current from a top-down (super-) 
gravity model, rather than in a phenomenological bottom- up AdS/QCD model. 
Moreover, the degrees of freedom in the dual field theory are identified explicitly. 

1. FIELD THEORY INTERPRETATION 

In this letter we consider 3+1-dimensional Af = 2 supersymmetric Yang- 
Mills theory, consisting of a Af = 4 gauge multiplet as well as two Af = 2 
supersymmetric hypermultiplets, called (0„, ip u ) and ((j)d,ipd)- If the masses are 
degenerate, the theory has a global U (2) flavor symmetry. 

A finite isospin chemical potential /i is introduced as the source of the oper- 
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ator 

Jq cc i>T 3 -f Q ip + 4>T 3 d (f> = n u - n d , (1) 

where n u / d is the charge density of the isospin fields, (4> u , 4>d) — 4> and (ip u , ipd) — 
ip. t 1 are the usual Pauli matrices. A non-zero vev ( Jq) introduces an isospin 
density as discussed in Q . The isospin chemical potential /i explicitly breaks the 
U{2) ~ U{1) B x St7(2)j flavor symmetry down to U{1) B X f7(l) 3 , where U(l) 3 
is generated by the unbroken generator r 3 of the SU(2)j. Under the U(l)s 
symmetry the fields with index u and d have positive and negative charge, 
respectively. 

In this letter we show that above a critical value for the isospin density the 
system is stabilized by a g-meson superfluid, i.e. a state with a non-vanishing 
vev of the flavor current component 

Jl cx-0t 1 73'0-|-^t 1 <93 ^ 
= 4>u 73 ipd + 4>d 73 ipu + bosons . 

The vev has p-wave symmetry and breaks both the SO(3) rotational symmetry 
as well as the remaining Abelian C7(l) 3 flavor symmetry. The rotational £0(3) 
is broken down to 50(2)3, which is generated by rotations around the x 3 axis. 
These symmetries are spontaneously broken. However only the J7(l)3, which 
acts on the gauge field on the D7 brane probe, is a dynamical symmetry in our 
approach, since we do not consider the backreaction of the D7 brane gauge field 
on the metric. This implies that we hold the metric background and thus also 
its symmetries fixed. Consequently only one Nambu-Goldstone boson is visible 
in our approach, namely the one due to spontaneous breaking of the dynamical 
U(l) 3 . 

Later we will compare our results to properties of superconductors: the 
condensate (J|) is the analog to the Cooper pairs of the BCS theory. Further, 
the £7(1)3 is the analog to the U(l) em and the current J 3 to the electric current 

Jem- 



2. HOLOGRAPHIC SETUP 

We consider asymptotically AdSs x S 5 spacetime which is holographically 
dual to N = 4 super- Yang-Mills theory with gauge group SU(N C ). The dual 
description of a finite temperature theory is an AdS black hole [H . We use the 
coordinates of [§], 

dfi2 = 2^ {-J dt2 + /df2 ) + f + ' (3) 

with 

/ = l-f • /=l + f ■ (4) 
where R is the AdS radius and dfl^ the metric of the unit 5-sphere. 
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Flavor degrees are added by embedding Nf coincident probe D7 branes with 
Nf <C N c into the AdS black hole background [13, EH, EH- Here we restrict 
ourselves to the case of zero quark mass and only comment on our massive 
results. In the massive case the physical interpretation does not change. For 
zero quark mass the induced metric G on the D7-branes is 

ds2(G) = ^ {-J df + /df2 ) + f ^ + RW * ■ (5) 

The Dirac-Born-Infeld (DBI) action determines the embedding of D7-branes as 
well as the gauge field on these branes, 

Sdbi = -T D7 Jd a £ Street (G + 2to'F)| , (6) 

with non-Abelian gauge field F = AA + cA -1 / 2 [A, A] on the D7-branes. The 
symmetrized trace prescription in this DBI action is only valid to fourth order 
in a'. However the corrections to the higher order terms are suppressed by 
NJ 1 . Here we use two different approaches to evaluate (HJ). First, we expand 
the DBI action to fourth order. Second, to make the calculation of the full 
action (|6|) feasible, we modify the symmetrized trace prescription by omitting 
the commutators of the generators t % and setting (t 1 ) 2 = 1. The numerics 
presented are obtained by the second prescription. However for both approaches 
we find properties of a superconductor, the critical exponent of one half up to an 
error of ten percent and a second order phase transition. The numerical values of 
thermodynamical quantities differ in the two approaches due to the modification 
of the effective action for the meson interaction. Still, the qualitative behavior 
remains the same. 

In the following we restrict to Nf = 2 flavors. An isospin chemical potential 
fi may be introduced by a non- vanishing time component of the non-Abelian 
gauge field A = Aq(q) t 3 on the D7-brane probe 0. The operator Jg (see ([2])) 
is dual to the non- vanishing gauge field component A\{q) on the D7 branes. 

The equations of motion for Aq(q) and A\(g) obtained from the DBI action 
([6]) are only satisfied if the asymptotic expansion near the boundary is of the 
form 



A 3-.._ _!§_££.... A 1 — + f7) 

2na' p 1 2-Kct' - 



where p = g/ Qh is the dimcnsionless AdS radial coordinate, qh is the radius of 
the horizon. According to the standard AdS/CFT correspondence, the absence 
of a constant term in the expression for A\ above ensures that the corresponding 
symmetry t/ (1)3 is broken spontaneously and not explicitly. Moreover, accord- 
ing to AdS/CFT, fi is the isospin chemical potential and the parameters d are 
related to the vev of the currents J by 

J3 = 2f(jg) & = 2f(4) 

° NfN c y/\T 3 ' 3 NfN c y/\T 3 ' 
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Figure 2: Dimensionless grand potential W7 (large plot) and specific heat C7 (small plot) 
versus temperature T/T c . Blue: the unstable phase with = 0. Red: the new stable phase 
with A3 ^ 0. Note that the dimensionful specific heat C7 oc T 3 Cj is rising above T c . 



2.1. Legendre transformation 

We Legendre-transform the action which gives rise to first order equations 
of motion which are easier to handle numerically. The conjugate momenta are 
defined by 

3 _ SSbbi 1 _ SSdbi , q n 

P °-5(d e Al)> P *-5{d e A\)- ™ 

These conjugate momenta are functions of the AdS radial coordinate g — and 
not constant in contrast to 0, 0, [HI — due to the non-Abelian term A|Aq in 
the DBI action. 

Using the asymptotics ([7]) , the densities d may be calculated from the asymp- 
totic value of the dimensionless conjugate momenta p — p / (2ira' N fTpj g^ ), 
Pq — > d$ and p\ — > — d\. The Legendre-transformed action is given by 



Sdbi — —NfTr>7 



k 2 p 4 P J V ' p 6 / 3 P 6 ff 2 

(10) 

The action simplifies to the result of [2| if A\ = 0. A non-zero gauge field A\, 
which also implies a non-zero conjugate momentum p\, decreases the value of 
the Legendre-transformed action S'dbi- Since the contribution of the flavors to 
the free energy F 7 is given by 

F 7 = TSdbi , (11) 

the phase with non-zero gauge field A\ is thermodynamically favored. Since the 
action S'dbi diverges on-shell due to the infinite AdS volume, the action must 
be renormalized. Details may be found in [14 |2I|. 



3. PHASE TRANSITION 

In this section we study the second order phase transition to the phase 
with non-zero gauge field A3. We work in the grand-canonical ensemble. The 
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contribution of the flavors to the grand potential is given by the on-shell action 
Sdbi 

n 7 = ™^w 7 = r Sl> B,, (12) 

where the dimensionless quantity WV is plotted in fig. [2l In this figure we also 
present the contribution of the flavors to the specific heat Cj , 

XN c N f V 3 T 3 d^ 7 
Cr- gg C ?-T— . (13) 

The temperature scale is defined by 

- 

T jl c I (do)c\ 3 



(14) 



with the dimensionless critical chemical potential Jx c w 2.85 and density (d(j)c ~ 
20.7. Fig.[5]shows a smooth transition in the grand potential and a discontinuous 
step in the specific heat at the phase transition. This implies that the transition 
is second order. 

From a fit to the numerical result for the order parameter d\ near the critical 
density, we obtain the critical exponent of the transition to be 1/2 up to an 
error of 10%, in agreement with the mean field theory result. Since the order 
parameter d\ and the density dp increase rapidly as T — > 0, we expect that the 
probe approximation ignoring the backreaction of the D branes to the geometry 
breaks down near T = 0. 

Let us compare our results to QCD. In QCD, the pion condensate is of 
course the natural state in isospin asymmetric matter. The condensation of a 
particle sets in if the isospin chemical potential is larger than the mass of this 
particle. According to this rule, the pions condense first in QCD since they 
are the Nambu-Goldstone bosons of the spontaneous chiral symmetry breaking 
and therefore the lightest particles. However the dual field theory which we 
consider in this letter is supersymmetric at zero temperature and therefore chiral 
symmetry cannot be broken spontaneously. In this supersymmetric theory, the 
vector and scalar mesons have the same mass at zero temperature. Due to 
finite temperature effects, the mass of the vector and scalar mesons can become 
different as we increase the temperature. It is a priori unclear which particle 
will condense. In our model we checked that the vector mesons condense first 
such that the p-meson condensation state, which we consider in this letter, is 
the physical ground state of our system near the phase transition. 

Since the system shares similarities with superconductors, it is interesting 
to investigate the remnant of the Meissner-Ochsenfeld effect. Switching on a 
magnetic field H 3 = Ff 2 along the lines of [HI, 0], we find a critical line in 
the magnetic field and temperature plane. This line separates the normal phase 
from the superconducting one at finite external magnetic field. Similarly to 
the discussion of [8], this suggests the existence of a Meissner phase at small 
magnetic field and temperature. 
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Figure 3: Sketch of our string setup: Strings spanned from the horizon of the AdS black hole 
to the D7-branes (green and blue plane) induce a charge at the horizon Q, Q, |17| . In the 
setup considered here, there are also D7-D7 strings present as shown in the figure. These D7- 
D7 strings are distributed along the AdS radial coordinate g, since they have to balance the 
flavorclectric and gravitational, i.e. tension forces. Thus these D7-D7 strings distribute the 
charges along the AdS radial coordinate, leading to a stable configuration of reduced energy. 
This corresponds to a superconducting condensate given by the Cooper pairs. 

4. STRING THEORY PICTURE 

In the string context, the non-zero fields Aq and A\ induce two non-zero 
flavorelectric SU{2) fields E\ = F ( % = A%A\ and E 3 g = F 3 g = -d e A% as well 
as a non-zero flavormagnetic field B\ g — F^ g = — d e A\ on the D7-branes (we 
use the notation of [9j, [4] with g being related to the radial AdS direction). 
These fields are generated by D7-D7 strings stretched between the two probe 
branes and by strings stretched from the D7-branes to the horizon. The D7-D7 
strings are new to our setup. They move from the horizon into the bulk and 
thus distribute the isospin charge. This stabilizes the system. 

Let us first describe the unstable configuration in absence of the field A\ . As 
known from [l7L Q , the non-zero field Aq is induced by fundamental strings 
which are stretched from the D7-brane to the horizon of the black hole. Since 
the tension of these strings would increase as they move to the boundary, they 
are localized at the horizon, i. e. the horizon is effectively charged. By increasing 
their density, the charge on the D7-brane at the horizon and therefore the energy 
of the system grows. In Q, the critical density was found beyond which this 
setup becomes unstable. In this case, the strings would prefer to move towards 
the boundary due to the repulsive force on their charged endpoints generated 
by the field E 3 g . 

The setup is now stabilized by the new non-zero field A\. This field is 
induced by D7-D7 strings moving in the x 3 direction. This may be interpreted 
as a current in x 3 direction which induces the magnetic field B\ = — d g A^. 

Let us now explain the mechanism by which the D7-D7 strings may prop- 
agate into the bulk. Due to the non-Abelian structure, the field E 3 and the 
magnetic field B\ Q induce the field i?| , corresponding to an interaction between 
the two string types. This field E% stretches the D7-D7 strings in the a; 3 di- 
rection. The position of the string in the g coordinate is fixed such that the 
gravitational force induced by the change in tension balances the force induced 
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Figure 4: Real part of the conductivity Re a versus frequency w = uj/(2ttT) at different 
temperatures: T w 0.90T C (black), T « 0.66T C (green), T w 0.46T C (blue), T fa 0.28T C (red). 



by the field Bf. This means that the energy of the setup is minimized. Our 
numerical calculations show that this is the case. 

The double importance of the D7-D7 strings is given by the fact that they are 
both responsible for stabilizing the new phase by lowering the charge density, as 
well as being the dual of the Cooper pairs since they break the t/(l)3 symmetry. 

5. FLUCTUATIONS 

The full gauge field A on the branes consists of the field A and fluctuations a, 

A = A 3 r 3 dt + A^dxs + a^T a dx" , (15) 

where r a are the SU(2) generators. The linearized equations of motion for the 
fluctuations a are obtained by expanding the DBI action in a to second order. 
We will analyze the fluctuations a 3 and X = a\ + ia| > Y = a\ — \a\. 

5.1. Fluctuations in a| 

We calculate the frequency-dependent conductivity a(u>) using the Kubo 
formula, 

a(u) = -G R (u 1 ,q = 0), (16) 

where G R is the retarded Green function of the current JS dual to the fluctua- 
tion a|> which we calculate using the method obtained in [18j. The current Jf is 
the analog to the electric current since it is charged under the U(l)s symmetry. 
In real space it is transverse to the condensate. Since this fluctuation is the 
only one which transforms as a vector under the SO (2) rotational symmetry, 
it decouples from the other fluctuations of the system. The real part of the 
frequency-dependent conductivity Re <j(uj) is presented in fig. UJ It shows the 
appearance and growth of a gap as we increase the condensate d\. Using the 
Kramers-Kronig relation, which connects the real and imaginary part of the 
complex conductivity, we find a delta peak at ui = in the real part of the con- 
ductivity, Re <y(uj) ~ irn s S(uj). As expected from Ginzburg-Landau theory, our 
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Figure 5: Movement of quasinormal modes at increasing temperature T: The different colors 
indicate the different fluctuations X, and o|. A Higgs mechanism is evident as explained 
in the text. 



numerics show that the superconductive density n s vanishes linearly at the crit- 
ical temperature, n s oc (1 — T/T c ) for T w T c . As a second distinct effect fig. 0] 
shows prominent peaks which can be interpreted as mesonic excitations, as con- 
firmed by our massive calculation. This is reminiscent of results for condensed 
matter systems where prominent quasiparticle peaks appear (e.g. [191]). 

5.2. Fluctuations in X = a\ + \a\. Y = a\— \a\ 

As shown in fig. 03 our setup is stable with respect to the fluctuations X 
and Y. Furthermore, fig. [5] shows that the quasinormal modes of higher ex- 
citations n > 1 move to larger frequencies and closer to the real axis. This 
corresponds to the formation of stable massive mesons. Such a behavior is 
known in gauge/gravity duality for mesons which are built from massive quarks 
(e.g. 0). Thus we observe a dynamical mass generation for the mesons which 
we expect to be similar to the Higgs mechanism in the field theory. This is dual 
to the gravity gauge fields eating the Nambu-Goldstone bosons. This reason- 
ing also applies to the massive mesonic excitations in the fluctuations a\. The 
explicit identification of the Nambu-Goldstone bosons and the massive vector 
mesons is postponed to future investigations. 

6. DISCUSSION AND OUTLOOK 

We found a stringy realization of holographic superconductivity for which 
the field theory action is known. It will be interesting to study the D7-D7 string 
condensate further, e. g. the drag force of strings pulled through this condensate. 
We expect that in the condensate this force vanishes in contrast to the original 
black hole background [2(J. The methods presented in this letter can also be 
applied to D2/D6 or D3/D5 (see also [HI]) systems. 
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